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Abstract 

We consider the bipartite matching model of customers and servers introduced by Caldentey, 
Kaplan, and Weiss (Adv. Appl. Probab., 2009). Customers and servers play symmetrical 
roles. There is a finite set C, resp. 5", of customer, resp. server, classes. Time is discrete 
, ^ and at each time step, one customer and one server arrive in the system according to a joint 

probability measure /x on C x S, independently of the past. Also, at each time step, pairs of 
matched customer and server, if they exist, depart from the system. Authorized matchings 
are given by a fixed bipartite graph {C,S,E C C x S). A matching policy is chosen, which 
decides how to match when there are several possibilities. Customers/servers that cannot be 
^-H matched are stored in a buffer. The evolution of the model can be described by a discrete 

time Markov chain. We study its stability under various admissible matching policies includ- 
I I ' ing: ML (Match the Longest), MS (Match the Shortest), FIFO (match the oldest), priorities. 

There exist natural necessary conditions for stability (independent of the matching policy) 
^> defining the maximal possible stability region. For some bipartite graphs, we prove that the 

. stability region is indeed maximal for any admissible matching policy. For the ML policy, we 

^ prove that the stability region is maximal for any bipartite graph. For the MS and priority 

I I policies, we exhibit a bipartite graph with a non-maximal stability region. 
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1 Introduction 



In queueing theory, customers and servers play different roles. Customers arrive in the system, 
accumulate in a buffer, get served by a server, and eventually depart. Servers on the other hand 
alternate between idle and busy periods but remain forever in the system. 

Within this framework, many variations and refinements are possible. For instance, we may 
consider a model with multi-class customers and distinguishable servers. A customer of a given 
class c must choose its server from a specified subset 5(c) of the servers. And of course, the subsets 
5(c) may intersect, see Figure [ij 




Figure 1: Queueing model of a call center. 

In this paper, we consider a model with the same multi-class flavor, but in which, by contrast, 
customers and servers play completely identical roles. We now argue that this simple symmetry 
requirement leads in a natural and ineluctable way to the bipartite matching model. 

By symmetry, both customers and servers should arrive into the system and depart from it. 
More specifically, upon completion of a service, both the customer and the server should depart 
simultaneously. To model arrivals, we have a priori more flexibility, but there is basically one 
non-trivial choice which is to assume that time is discrete and that customers and servers arrive 
in pairs. 

Consider indeed the simplest possible model with continuous-time arrivals: (i) there is only 
one class of customers and one class of servers; (ii) customers, resp. servers, arrive according to 
a Poisson process of rate A, resp. /i; (iii) services have duration 0. Let us describe the state by 
Z = X — Y , where X is the number of unmatched customers and Y the number of unmatched 
servers. The process Z is a birth-and-death continuous-time Markov process on Z with drift A — /i. 
It is either transient (if A 7^ /i) or null recurrent (if A = /i), but it is never positive recurrent. 

Let us switch to discrete-time i.i.d. arrivals. At each time step, a batch of customers and a 
batch of servers arrive into the system. If the size of the batches are allowed to be different for 
customers and servers, then we are back to the continuous-time situation, and even the simplest 
model is never positive recurrent. 

Therefore to get a non-trivial model, the natural assumption is that exactly one customer and 
one server arrive into the system at each time step. The resulting model is symmetric in another 
respect: both arrivals and departures occur in pairs. 

For simplicity, we always assume that the service durations are null. So the model is specified 
by: (i) the finite set C of customer classes and the finite set 5 of server classes; (ii) the probability 
law /i on C X 5 for the arrivals in pairs; (iii) the bipartite graph (C, S, E C C x S) giving the 
possible matchings between customers and servers (hence the possible departures in pairs); (iv) 
the matching policy to decide how to match when several choices are possible. We consider so 
called admissible policies which depend only on the current state of the system. Under these 
assumptions, the buffer content evolves as a discrete-time Markov chain. 

We call this model the bipartite matching model. 

The bipartite matching model has been introduced by Caldentey, Kaplan, and Weiss [5] , under 
an additional assumption of independence between arriving customers and servers (Vc, s, /i(c, s) — 
n{c, S) ^{C , s)) , and for the FIFO policy. In their paper, the authors mention several possible 
domains of applications ranging from call centers to crossbar data switches. They also provide 
references to papers on related models. We refer the interested reader to jSj for details. 
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In the bipartite matching model, there is an equal number of customers and servers at any 
time. But the matching constraints may result in instability with unmatched customers and servers 
accumulating. It turns out that proving stability, i.e. positive recurrence of the Markov chain, is 
highly non-trivial. Given a bipartite graph {C,S,E), there exist natural necessary conditions on 
/i for stability to hold true. When these conditions are also sufficient, we say that the stability 
region is maximal. 

Caldentey & al conjecture that any bipartite graph has a maximal stability region for the 
FIFO policy [H Conjecture 4.2]. They prove the conjecture for some specific models (under the 
additional assumption of independence of customers and servers): (i) the N model defined by 
C = {1, 2}, S = {!', 2'}, E = C X S - {{2, 2')}, (ii) the W model, i.e. the matching model version 
of Figure [I] (iii) the NN model of Figure [2] In the first two cases, they are also able to compute 
explicitly the stationary distribution. For the last case, the proof of stability is already intricate. 

In the present paper, we consider the stability issue for various admissible matching policies: 
ML (Match the Longest), MS (Match the Shortest), FIFO (match the oldest), random (match 
uniformly) , priorities. The irreducibility of the Markov chain describing the model is not granted, 
and we first study this question in detail (Section |4]). Then we obtain the following results: 

- sufficient conditions under which any admissible policy is stable (Section [5|; 

- for the NN model, the MS policy and some priority policies do not have a maximal stability 
region (Section |6]); 

- for any bipartite graph, the ML policy has a maximal stability region (Section [?]). 

We do not know if the stability region is always maximal for the FIFO and random policies. 

Notations. Denote by N = {0, 1,2,.. .} the set of non-negative integers. Let A* be the free 
monoid generated by A. For any word w e A* and any B C A, set \w\b — #{i | Wi £ B}, the 
number of occurrences in w of letters from B. For B = {b}, we shorten the notation to \w\b. 
Furthermore, for any w G A* , set [w] := (|w|a)aeA (the commutative image of w). 

2 The bipartite matching model 

We now proceed to a more formal definition of the model. 

Definition 2.1. A bipartite matching structure is a quadruple {C,S,E,F) where 

• C is the non-empty and finite set of customer types; 

• S is the non-empty and finite set of server types; 

• E d C X S is the set of possible matchings; 

• F d C X S is the set of possible arrivals. 

The bipartite graph (C, S, E) is called the matching graph. It is assumed to be connected. The 
bipartite graph (C, S, F) is called the arrival graph. It is assumed to have no isolated vertices. 

The two assumptions in Def . |2.1| are made without loss of generality, see Remarks [l] and [2j 

In Figure [2] we give an example of a matching graph with 3 customer and 3 server types, called 
the "NN graph" in the following. 




2' r 

Figure 2: NN graph. 
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Customers and servers play symmetrical roles in the model. Also E and F play dual roles. 
The graph (C, S, E) defines the pairs that may depart from the system, while the graph (C, S, F) 
defines the pairs that may arrive into the system. 

Definition 2.2. A bipartite matching model is a triple [(C, 5, i?, i^), /x, Pol], where 

• {C,S,E,F) is a bipartite matching structure; 

• fi is a probability measure on C x S satisfying 

supp{fi) = F, suppi^ic) = C, supp{ns) = S , (1) 

where o,''^d /ig are the C and S marginals of fi. 

• Pol is an admissible matching policy (to be defined in ^2.S^ . 

Observe that we can simplify the notation to [{C, S, E), fj,, Pol]. We say that the model 
[(C, S, E),iJ,, Pol] is associated with the structure (C, S, E, F). 

Remark 1. For (fl) to have solutions, (C, S, F) must be without isolated vertices, the assumption 
made in Definition 2.1| This is not a real restriction: if it is not satisfied, we can consider a new 
model without such customer or server classes. 

A realization of the model is as follows. Consider an i.i.d. sequence of random variables of law 
/i, representing the arrival stream of pairs of customer/server. A state of the buffer consists of 
an equal number of customers and servers with no possible matchings between the classes. Upon 
arrival of a new ordered pair (c, s), two situations may occur: if neither c nor s match with the 
servers/customers already present in the buffer, then c and s are simply added to the buffer; if 
c, resp. s, can be matched then it departs the buffer with its match. If several matchings are 
possible for c, resp. s, then it is the role of the matching policy to select one. An admissible policy 
selects according to the current state of the buffer (and not according to the whole history of the 
buffer contents, for instance). The resulting evolution of the buffer is described by a discrete-time 
Markov chain. 



2.1 State space description 

Depending on the matching policy, we consider either a commutative (e.g. for Random) or a 
non-commutative (e.g. for FIFO) state space description. The different policies considered in the 
paper will be formally defined in §2.2[ 

Let us choose a matching graph {C,S,E). We introduce the following convenient notations: 
C(s) is the set of customer classes that can be matched with an s-server; S{c) is the set of server 
classes that can be matched with a c-customer: 

S{c) ^{seS : (c, s) e E}, C(s) = {c G C : (c, s) e E}. 

For any subsets A C C, and B C S, we define 

S{A) = UceAS{c), C{B) = \JseBC{s). 



Commutative state space. A state of the system is given by {x,y), x — {xc)c£C and y — 
iys)seSj where Xc denotes the number of customers of type c and y^ the number of servers of type 
s. The commutative state space is: 

£-{(.x,t/)eN^xN^ : = V(c,s)eS, = 0}. (2) 

Non-commutative state space. A state of the system is given by two finite words of the same 
size fc > 0, respectively on the alphabets C and S, describing unmatched customers and servers. 
The non- commutative state space is: 

£^[{u,v)eUk>o{C''-^S'') : (M,H) belongs to (Ell. (3) 
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Facet. Both the commutative and the non-commutative state space can be decomposed into 
facets, defined only by the non-zero classes. 



Definition 2.3. A facet is an ordered pair {U, V) such that: U C C,V C S and U x V C 

(C X S — E). The zero-facet is the facet (0, 0), we denote it shortly by 0. 

For a facet T = {U,V), define: 

C@{T) = C{V), S@{T) = S{U), 

CoiJ") = C- (C.(J-) U C®(J-)), So{F) = S- {S,{F) U S@{T)). 

We alleviate the notations to C,,S,, C^, . . . , when there is no possible confusion. The symbol 
• stands for the non-zero classes, the symbol ® for the classes that are forced to be at zero (since 
they are matched with non-zero classes), and the symbol o for the classes that happen to be at 
zero. 

The following notion will play an important role later on. 
Definition 2.4. A facet T is called saturated if Co {J') = % or 5o(J^) = 0. 

In Figure [3) the facet on the left is non-saturated, while the one on the right is saturated. 



123 123 




Figure 3: NN graph: facets ({3}, {3'}) and ({2}, {3'}). 



Graphical convention. A facet F can be represented graphically by coloring the nodes of the 
bipartite graph according to the above convention (see Figure [3] for an illustration): 

- nodes in C,{J-) and S,{J-) are represented as filled circles; 

- nodes in C@(J^) and S@{J^) are represented as double circles; 

- nodes in Co{J^) and So{J^) are represented as simple circles. 

In Figure |4] we have represented the facets of the NN graph. The more complex case of the 
NNN graph will be given in Section [5] Figure [13] 








{1},{1') {2),{3'} {3},{3'} {3},{2'} 




{2,3}, {3'} {3},{2',3'} 
Figure 4: Facets for the NN graph. 

Algorithm [1] takes as input a matching graph and returns as output the set of facets. The 
termination and correctness of the algorithm are easily proved. 
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Algorithm 1: Comijutalioii of tin" facets 



Data: A bipartite graph G = {C, S, E). 
Result: Facets - set of all facets of G. 
begin 

Facets ^ 0; New ^ 0; 

foreach G G x S - E do New ^ New U {{{i}, {j})}; 
while New ^ do 

Facets -s— Facets U New; 
Old ^ New; New ^ 0; 
forall the H,K. G Old such that H ^ JC do 
if C,{n) = C.(/C) or S,{n) = S'.(/C) then 
Z ^ (C.CH) U C.(/C), S,{n) U S,{IC)); 
New <(- New U {Z}; 

Facets ^ Facets U {0}; 
return Facets; 



2.2 Admissible matching policies 

Informally, a matching policy is admissible if: 

- only the current state of the buffer is taken into account; 

- priority is given to customers/servers that are already present in the buffer: if the state is 
{u, v) and the new arrival is (c, s) G E, then c and s are matched together iff there are no 
servers from S{c) in v and no customers from C(s) in u. 

It results from the first point that an admissible matching policy can be described as a mapping 
Q : £ X {C X S) ^ £ which returns the new state of the system after an arrival. The second point 
is called the buffer-first assumption. It is not a real restriction: a matching policy that always 
gives priority to new arrivals can be seen as a special case of the above with an arrival probability 
/X such that n{E) = 0. 

We now define admissible policies formally, distinguishing between the non-commutative and 
commutative state spaces. 

For a word w G A'^ and i G {1, . . . , k}, we denote by := wi. . . Wj-iWj+i . . .Wk the subword 
of w obtained by deleting Wi. 

Definition 2.5 (Non-commutative case). A matching policy is admissible if there are functions 
$ and \1/ such that: 



{uc,vs), 


/ Hc{s) 


= 0, 


kls(c) 


= 0, (c. 


s)^E 




f W\c{s) 


= 0, 


\v\s(c) 


= 0, (c. 


s)gE 


(U[*(ti,s)],^^[*(t;,c)])) ' 


f Hc{s) 




kls(c) 


y^O 






f Hc{s) 


^0, 


kls(c) 


= 






f VAc(s) 


= 0, 


l«|s(c) 


^0 





The FIFO and LIFO policies are admissible matching policies with functions $ and ^ as 
follows: 

• FIFO : s) = argmin{?ife e C'(s)}, ^{v,c) = argminjwfc e S{c)}. 

• LIFO : <I>(m, s) — argmaxjufc e G{s)}, 5'(w,c) = argmaxjufc G S{c)}. 

For c G C, let Cc G N*-^ be defined by (ec)c = 1 and (ec)d = 0,d ^ c. For s G S, let be 
defined accordingly. 
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Definition 2.6 (Commutative case). A matching policy is admissible if there are functions <& and 
^! such that: 



{x, y) (c, s) = < 



{x + ec,y + Cs), if Xc(s) = 0, ys{c) = 0, (c, s) E 

{x,y), if xc{s) = 0, ys{c) = 0, (c,s) e E 

{x - e^i^^^s),y - e^{y,c)), if xc(s)i^^i Vsic) ^ 

{x - e^i^^^s) + Gc, y), if xc{s) 7^ 0, ys{c) = 

, {x,y-e^i^y^c)+es), if xc(s) = 0, ys(c} ^ 



The following commutative matching policies are admissible (for RANDOM, ML, and MS 
policies $(m, s) and ^'(i;,c) are random variables): 

• PR (Priorities). For each customer type c e C, we define a priority function : S{c) — )■ 
{!,..., |S'(c)|}. Similarly, for each server type s € S", we define Ps ■ C'(s) — > {!,..., |C(s)|}. 
In the case of several matching options, a customer/server is matched with the server/customer 
that has the highest priority (greatest value of the priority function) . It is convenient to spec- 
ify the priorities by two |C| x |S'| matrices A and B defined by: 

ads), {c,s)eE ^ _//3,(c), {c,s)eE 



0, otherwise 1 0, otherwise 

Then $(a;,s) ~ arg max{/3s (c) : c G C(s),Xc > 0} and '^{y,c) — argmax{ac(s) : s € 
5(c), >0}. 

• RANDOM : ^{x,s), resp. ^'(?/,c), is a random variable valued in C(s), resp. S{c), and 
distributed as ( X^iecfs) ) ' resp. ( ipsfc) ) ■ Intuitively, the match 
is chosen uniformly among all possible ones. 

• ML : ^(x,s), resp. ^(j/,c), is a random variable uniformly distributed on argmax{a:i : i € 
C(s)}, resp. argmax{yi : i G S{c)}. 

• MS : $(x,s), resp. ^'(?/,c), is a random variable uniformly distributed on argminjxi > : 
i £ C(s)}, resp. argmin{?/,i > : i £ S{c)}. 



3 Necessary conditions for stability 

To introduce the main ideas, consider first a simpler finite and deterministic problem. Let (C, S, E) 
be a matching graph. Consider a batch of customers x E N"^ and a batch of servers y E of 
equal size: J^c = J2s Vs- ^ perfect matching of x and y is a tuple m E such that: 

Vc E C, Xc = ^ nics, Vs E S, ys = ^ nics ■ 
ses{c) cec{s) 

By Hall's Theorem (aka the "marriage Theorem"), there exists a perfect matching if and only if: 

Ecec/ < J2seS{U) Vs^ VC/ C C 
J^sevVs <J2ceciv)^c, VV c S 

A perfect matching, if there is one, can be obtained by restating the model as a flow network and 
by solving the maximum flow problem for which efhcient algorithms exist [H |3] . 

The bipartite matching model is much more complicated: first it is random, and second the 
matchings have to be performed on the fly, at each time step. However the two ingredients of 
the simpler model will play an instrumental role in the analysis: (i) the conditions NCOND, to 
be defined in ([s]), are related to Q; (ii) the restatement as a flow problem is used in most of the 
proofs. 
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Consider now a bipartite matching model [(C, S, E),fi, Pol]. Wc identify the model with the 
Markov chain on the state space £ describing the evolution of the buffer content. 

Let P be the transition matrix of the Markov chain. A probability measure tt on £ is stationary 
if ttP = TT. It is attractive if for any probability measure v on the sequence of Cesaro averages 
of uP^ converges weakly to tt. 

Definition 3.1. The model is said to be stable if the Markov chain has a unique and attractive 
stationary probability measure. 

It implies in particular that the graph of the Markov chain has a unique terminal strongly 
connected component with all states leading to it. 

Let fic be a probability measure on C and ^J-s ^ probability measure on S. Define the following 
conditions on (ficTfJ-s)' 



NCOND 



^isiv) < ^iciciv)), yvcs 



(5) 



The above conditions appear in [2]. They have a natural interpretation. Let nc and /is be the 
marginals of the arrival probability /i. Customers from U need to be matched with servers from 
S{U). The first line in NCOND asks for strictly more servers in average from S{U) than customers 
from U. The second line has a dual interpretation. 

Using the Strong Law of Large Numbers, we also see that the arrivals up to time n satisfy Q 
for all values of n large enough if and only if NCOND is satisfied. 

Lemma 3.2. The conditions NCond are necessary stability conditions: if the Markov chain is 
stable then the conditions NCOND are satisfied by the marginals of fi. 

Proof. We suppose that the conditions NCond are not satisfied. 

Assume first that there exists U C C such that fJ,c{U) > /is(S'([/)). Let An, resp. i3„, be 
the total numbers of customers of type U, resp. servers of type S(U), to arrive in the system up 
to time n. Let A"„ be the number of customers of type U present in the system at time n. By 
definition, Xn > An — Bn. By the Strong Law of Large Numbers, we have, a.s.. 



lim ■ 



Ar, 



Hc{U), lim — 
n n 



So the Markov chain is transient. Similarly, if there exists V C S such that usiV) > ^c{C{V)), 
the model is unstable. (This part of the argument appears in ^ Prop. 3.4].) 
Assume now that there exists U C C,U ^ C, such that 



(6) 



Observe that S{U) ^ S, otherwise we would have Hc{U) — IJ-s{S) — 1 which would contradict 
U =^C. Set y = S-S{U). Eqn ^ is equivalent to: ns{V) ^ fic{C-U). The bipartite matching 
graph (C, S, E) is represented in Figure [5] By assumption, U 'k V D E = %. 



C 


u 


c-u 










s 


S(U) 


V 



Figure 5: The bipartite graph (C, S, E). 



We have 



^l{Uy.V) - iic{U) - ^l{U y. S{U)) 
fi{{C-U)x S{U)) = fis{S{U))^fi{UxS{U)) 
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Using Eqn we get 

m({7 xV)^ ^l{{C -U)x S{U)) . (7) 

Let Dn be the number of departures of type {C — U) x S{U) up to time n. Let A„, i?„, and 
Xn be defined as above. Set Z„ = An — Bn- For an arrival of type U xV , the Z-process makes 
a +1 jump, for an arrival of type {C ~ U) x S{U), the Z-process makes a -1 jump, otherwise the 
Z-process remains constant. We have X„ > — {Bn — Dn) > Z^- We have two cases: 

- If pl(\J X F) > 0, then, according to ([t]), the Z-process is null recurrent. 

- If X F) = 0, then for any initial condition such that > 0, a.s. Xn > Xq, Vn. 
Hence, in both cases, the model cannot be stable. □ 

Remark 2. Consider a non-connected matching graph {C,S,E). Consider a probability ^ and 
an admissible matching policy such that the bipartite matching model is stable. Let {C , S' , E') 
be a connected subgraph of (C, S, E). Following the exact same steps as in the proof of Lemma 



3.2 we prove that 

^Jic{C')^ lists'), ^i{C' X (S - S')) = 0, fi{{C-C')x S') = 

(otherwise the Markov chain is either transient or null recurrent). Therefore, we can decom- 
pose the model into connected components and treat them separately. Hence the assumption of 



connectedness of (C, S, E) in Def. 2.1 was made without loss of generality. 



3.1 Complexity of verifying NCOND 

Let us fix (C, S, E) and the probability measures (/xci Ms) such that supp(/i(7) — C, supp(^5) — S. 
We want an efficient algorithm to decide if the conditions NCOND are satisfied. 

The number of inequalities in NCond is exponential in |C| + IS*]. So checking directly if all 
the inequalities are satisfied is a method whose time complexity is exponential in |C| + \S\. To go 
beyond, we need additional material. 

We use the standard terminology of network flow theory, see for instance |4j. Consider the 
directed graph 

= (C U 5 U {z, /}, £; U {{i, c),c eC}U {{s, f),s e S}) . (8) 

Endow the arcs of E with infinite capacity, an arc of type (i,c) with capacity fic{c), and an arc 
of type (s,/) with capacity /is(s). 




Figure 6: The graph associated with the NN model of Figure [2j 

Recall that a cut is a subset of the arcs whose removal disconnects i and /. The capacity of a 
cut is the sum of the capacities of the arcs. Set A = EU {(z, c), c G C} U {(s, /), s G S}. Recall 
that T : A is a. flow if: (i) Vc, T{i, c) = EseS(c) T{c, s), Vs, Ecec(.) T{c, s) = r(s, /); (n) 

y{x,y) e E, T{x,y) is less or equal to the capacity of {x,y). The value of T is X]c^(*''^) ~ 

Let NCOND< be the set of inequalities obtained from NCond by replacing the strict inequal- 
ities by large inequalities. 
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Lemma 3.3. There exists a flow of value 1 in Af ijf [fJ-Ci fJ-s) satisfies NCond<. There exists a 
flow T of value 1 such that T(c, s) > for all (c, s) ^ E iff (/ic: Ms) satisfies NCOND. 



The first part of Lemma 3.3 is proved in [2, Prop. 3.7]. We repeat tlie argument for complete- 



Proof. The celebrated Max-flow Min-cut Theorem 4j states that the maximal value of a flow is 
equal to the minimal capacity of a cut. Observe that the set of arcs {(z, c), c e C} forms a cut of 
capacity 1. Therefore the maximal flow is < 1 and it is 1 iff all cuts have a capacity > 1. 

To be of finite capacity, a cut must contain only customer arcs {(i, c), c G C} and server arcs 
{(s, /), s e S}. Consider a subset A = {{i, c), c e Ci} U {(s, /), s e Si}. Set C2 = C - Ci and 
5-2 = 5'- Si. The set ^ is a cut iff C2 x 52 n £; = 0, equivalently iff 5(C2) C Si and C(52) C Ci. 
Also the capacity of A is ficiCi) + usiSi). 



ci 



C2 




Figure 7: Illustration of the proof of Lemma |3.3[ 

Assume that the cut {(z,c),c e Ci} U {{s,f),s e Si} is of capacity strictly less than 1. We 
have 

Aic(Ci) + ^Is{Sl) < 1 ^ ticiCi) < ^^s{S2) ■ 

But C(52) C Ci so, if NCOND< is satisfied, we must have: 

I^S{S2) < t^c{C{S2)) < ficiCl) . 

So we have proved that NCOND< is not satisfied. 

The other way round, if NCond< is not satisfied, then there exist Ci,52, C(52) = Ci such 
that Mc(Ci) < HsiS2). Set C2 = C- Ci and 5i = 5 - 52. By definition, C2 x S2 n E ^ H), 
therefore {{i, c), c G Ci} U {(s, /), s € Si} is a cut. Its capacity is ^c{Gi) + fJ-siSi) < 1. 

By contrapposing the above, we get that: 



NCOND< satisfied 



all cuts have a capacity > 1 



maximal fiow is 1 



We now prove the second part of the lemma. Assume that the conditions NCOND are not 
satisfied. If the conditions NCond< are not satisfied either, then by the first part of the proof 
there exists no flow of value 1. Assume now that the conditions NCOND< are satisfied. Then 
there exists U C C, C/ 7^ C, such that ^c{U) — iis{S{U)). Let T be any flow of value 1. Using 
the flow relation for L/, we get: 



{c,s)<^Uy.S(U) 



E T{^,c)=^lc{U) 

(i,c)g{i}x;7 



Using fJ.c{U) — ns{S{U)) and the flow relation for S{U), we deduce that: 



J2 T(c,s) = Ms(5(t/)) 
{c,s)euxsiu) 



ic,s)e{c-u)xS{u) 
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Now it follows from the connectedness of (C, S, E) that [C - U) x S{U) n E ^ ^. We conclude 
that the flow T is such that T(c, s) = for some (c, s) e E. 

Assume now that the conditions NCond are satisfied. Fix such that < < Consider 
the function : A — >• ]R_|_ defined by 

rj for {x,y) — (c, s) e E 

T^^ix^y) = <\S{c)\ 7] for = (i,c) 



\C{s)\7j for = 



By construction is a flow. Set 



~ / ^ /^c(c) - \Sic)\r] ^ fisjs) - |C(s)h , , 

For 7/ small enough, observe that JIq^ resp. /Is, is a probability measure on C, resp. on 5*. Choose 
T] small enough such that {Jlc, /^s) satisfies NCond. This is possible since the conditions NCond 
are open conditions. 

Consider the directed graph Af, see ([s]) , with new capacities on the customer and server arcs 
defined by Jlc and ptg. By applying the first part of the proof, there exists a flow T : A — > M+ of 
value 1. Define 

T:A^R+, T ^ T,, + {1 - \E\ri)T . 

By construction T is a flow for the graph M with the original capacity constraints (/ic for the 
customer arcs and fis for the server arcs). The value of T is 1 and it satisfies T{x,y) > for all 
{x, y) G E. This completes the proof. □ 

There exist algorithms to find the maximal flow which are polynomial in the size of the un- 
derlying graph, independent of the arc capacities. For instance, the classical "augmenting path 
algorithm" of Edmonds & Karp [3] operates in 0((|C| + |S'|)|£'p) time, and there exist more 
sophisticated algorithms operating in 0((|C| + \S\y^) time. 

Take one of these polynomial algorithms, call it MaxFlow and consider it as a blackbox. We 
build on this to design a polynomial algorithm to check NCond. Let us detail the construction. 

Lemma 3.4. Define {]lc,Jj-s) as in The pair (/ic/ig) satisfies NCond iff the pair {Jlc Jig) 
satisfies NCOND for rj strictly positive and small enough. 

Proof. Assume that (pcil^s) satisfies NCond. Since we are dealing with open conditions, any 
small enough perturbation of {iJ,c,t^s) still satisfies NCond. 

Assume now that {^CtI^-s) does not satisfy NCond. There exists U C C,U ^ C, such that 
t^ciU) > ^^siS{U)). By using we get 

{l-\E\r^)jlc{U)+(j2\Sic)\)v > {l-\E\rj)MSm + {T. V 

ceu ses{u) 

{1^\E\tj) Jlc{U) + \En{UxS{U))\ r, > {l ~ \E\r^) Jis{S{U)) + \E D {C{S{U)) x S{U))\ v ■ 

By definition, we have U C C{S{U)). We conclude that Jlc{U) > lis{S{U)). So the pair (JlcJis) 
does not satisfy NCOND. □ 



Using Lemmas 3.3 and 3.4 NCond is satisfied iff MaxFlow(A/', /ic, Ms) returns 1 for 77 small 
enough. So the trick is to run MaxFlow on the input (A/", mc, Ms) by considering 77 as a formal 
parameter made "as small as needed" . 

The precise meaning is the following. If xi, X2,yi,y2 G Hi, then: (xi + yirf) + {x2 + y2V) = 
{xi + X2) + {vi + y2)v- Furthermore, 

[xi + yiTj = X2 + t/2??] ■^=^ [xi = X2, yi = 7/2] 

[xi + yiT] <X2 + y2v] ■^=^ {xi < X2) or {xi = X2, yi < ^2) (10) 
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So 7] is small enough not to reverse any strict inequality. When running MaxFlow on (A/", Jlc, Jj-s), 
the algorithm deals with values of the type (x + yr/) , and adds and compare them according to the 
above rules. Now observe that the algorithm stops in finite time, so it will have performed only 
a finite number of operations. Therefore, it would be possible, a posteriori, to assign to 77 a value 
which would be small enough to enforce (10 1. 



Algorithm 2: Checking the necessary stability conditions 

Data: {C,S,E), {^C:IJ'S) such that suppdic) = C, supp{^s) — S. 

Result: "Yes" if NCOND, "No" if ^(NCond) 

begin 

Compute Af, Jlc, ]ls\ 

if MaxFlow(A/', ptcMs) = 1 then 
I Result ^ Yes; 

else 

|_ Result ^ No; 
return Result; 



The termination is obvious and the correctness follows from Lemmas 13.31 and 13.41 

Proposition 3.5. Given a bipartite model [[C,S,E)^^], there exists an algorithm of time com- 
plexity 0{{\C\ + 151)^) to decide i/NCOND is satisfied. 



4 Connectivity properties of the Markov chain 

Define the following property for the transition graph of the Markov chain: 

UTC : a unique (terminal) strictly connected component with all states leading to it. 



Property UTC is necessary for stability as defined in Def. 3.1 But property UTC is not granted 
in bipartite matching models and counterexamples are given below (Examples [2] and |3]) . In fact, 
we will see that we are in an unusual situation: the necessary stability conditions NCOND turn 



out to be sufficient conditions for the property UTC (Theorem 4.3 1 ! Observe also that property 
UTC is weaker than irreducibility, and we will give an example of a model satisfying NCOND and 
UTC without being irreducible (Example |4]). 

4.1 Stable structures 

To establish property UTC, we make a detour by introducing and studying a notion of independent 
interest: stable structures. 

Definition 4.1. A bipartite matching structure {C,S,E,F) is stable if there exists a probability 
measure /i satisfying |Ip and whose marginals fic o,^d fis satisfy NCOND. 

The justification for this terminology will appear in ^ we prove there that under the ML 
policy, any model satisfying NCOND is stable. So a structure is stable iff there exists an associated 
model which is stable. 

First of all, there exist stable structures. 

Example 1. Consider {C\S,E,C x S), where {C,S,E) is the NN bipartite graph of Figure [2j 

Let 

: Mc(l) = /^c(2) = 2/5, y.c{i) = 1/5, lis ■ Ms(l') = Ms(2') = 2/5, /is(3') = 1/5 . 

The product measure /i = x Ms has marginals /ic and and we check that (^CjMs) satisfy 
NCOND. Also it is easily proved that for any admissible matching policy, the graph of the Markov 
chain is irreducible. 
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On the other hand, there exist unstable structures. We illustrate this on two examples. 

Example 2. Consider the structure (C, S, E, F) where (C, S, E) is the NN graph of Figure[2j and 
where 

{(1,3'), (2, 2'), (3,1')}. 

Consider any /i with supp(/x) = F. We have /ic(l) — Ms(3') = /i(l,3') which violates NCond for 
V — {3'}. We can also prove that the property UTC is not satisfied. Consider a state of the type 
(x, y) with X = y = (0, 0, fc), for some fc > 0. Any one of the three possible arrivals leave the state 
unchanged. In particular, there is an infinite number of terminal components. 

Example 3. Consider the bipartite matching structure defined in Figure |8] The graph (C, S, E) 
is represented on the left of the figure, while the graph (C, S, F) is represented on the right. 




Figure 8: The matching graph (C, 5, E) on the left, and the arrival graph (C, F) on the right. 

Consider any /i with supp(/i) = F. We have 

fisUr, 2'}) = /i(3, 1') + fi{4, 2') < /ic({3, 4}) , 

which contradicts NCond for U — {3,4}. We can also prove that the property UTC is not 
satisfied. Consider a state {x,y) with ^3 + 2:4 = fc > 0. Reducing the number of customers of 
types 3/4 would require an arrival of type (1, 1') or (1, 2') or (2, 1') or (2, 2'). But none of these 
pairs belong to F. Therefore it is impossible to reach a state {x' , y') with ccg + 0:4 < X3 + X4. On 
the other hand an arrival of type (3,3') or (3,4') or (4,3') or (4,4') strictly increases the number 
of customers of types 3/4. Hence all the states are transient, and there is no terminal strongly 
connected component. 

Stability of a structure is a decidable property. There exists a probability measure /i with the 
requested properties iff the following system of linear inequalities in the indeterminates /^(c, s), c € 
C, s G 5, have a solution: 

X](c,s)GCxS ^('^' ^) ~ 

m(c,s) > 0, 
^ m(c,s) = 0, 
Mc(c) = EsesA*(c>5), 

NCOND . 

However, the number of inequalities is exponential in \C\ + We are going to propose a criterion 
which is much simpler, both conceptually and algorithmically. 

Consider a bipartite matching structure (C, 5, E, F). Define F — {{s,c) \ (c, s) € F}. Associate 
with the structure the directed graph (C U S', i? U i^), in other words the nodes are C U S" and 
the arcs are 

c — > s, if (c, s) E, s — > c, if (c, s) G F . 
We have represented in Figure |9] the directed graph associated with the structure of Example [3) 



V(c,s)eF, 

V(c, s)£CxS~F, 

Vs e S, 
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Figure 9: The directed graph associated with the structure of Figure [8j 

The graph of Figure [9] is not strongly connected: the four nodes on the right form a strongly 
connected component. Similarly, the directed graph associated with the structure of Example [2] 
is not strongly connected. On the other hand, the directed graph associated with the structure of 
Example [T] is strongly connected. This is not a coincidence. 

Theorem 4.2. Let (C, S, E, F) be a bipartite matching structure. The following two properties 
are equivalent: 

1. {C,S,E,F) is a stable structure; 

2. {C U S, E (J F) is strongly connected. 

In particular, one can decide if a structure is stable with an algorithm of time complexity 
Oi\C\\S\)hy testing the strong connectivity oi {C U S, E \J F). 

Proof of Theorem \4.^ Assume that {C,S,E,F) is a stable structure. Let ^ be a probability 
measure satisfying (|lj) and NCOND. Suppose that there exist c € C, s € S, with no directed path 
from c to s in (C U 5, -E U F). Let succ(c) be the set of nodes that can be reached starting from c 
in {CUS,EUF). Set 

Ci = C n succ(c), Si^Sn succ(c), C2 = C-Ci, S2^S-Si. 

By assumption, s G 5*2. The following two properties hold: 

Ai(C2,5i) = 0, {CixS2)nE = (l}. 

Using Ai(C2, ^i) = 0, we get 

^ls{Sl)^^l{Cl,Sl)<^lc{Cl). 

But using [(Ci X 52) n S = 0] and NCond for [/ = Ci, we get 

From this contradiction, we deduce that for all c e C, s G S, there exists a directed path from c 
to s in {C U S,E U F). Similarly, we can prove that for all s E S,c £ C, there exists a directed 
path from s to c in {C U S, E (J F). 

Assume now that (C U S,EU F) is strongly connected. Consider the matrices A G M|jl^'^ and 
B e Rf"^ defined by 

^^^^|l/l^(c)| ii{c,s)eE ^^^^ fi/#{d,(s,d)ei?} if(s,c)eF 

1 otherwise ' 1 otherwise 

Consider the matrix AB e M'j;^'^. By construction, we have {AB)cd > if and only if there is 
a path of length 2 from c to d in the graph {C U S,EU F). Since {C S.EVJF) is strongly 
connected, we deduce that AB is irreducible. Clearly the spectral radius of AB is 1. Applying 
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the Perron-Frobenius Theorem [F, we obtain the existence of a hne vector x £ M'j: such that: 
Vc, Xc > 0, X^c^c — li and xAB = x. Set y = xA. Define the probabiHty measure /i on C x by 
/x(c, s) = UsBsc- I3y construction, we have /Ltc = x, /xg = y. Also, by construction, supp(/i) = F. 
Define the function T : A ^ R+ hy 

Vc e C, T(i, c) = Xc, Vs G S*, T(s, /) = y^, V(c, s) e E, T{c, s) = x^^cs ■ 



By construction, T is a flow of value 1 such that T(c, s) > for all (c, s) G i?. Using Lemma 3.3 



we get that (x,y) = (mciMs) satisfies NCond. □ 



4.2 Back to property UTC 

We now have all the ingredients needed to prove the following result. 

Theorem 4.3. Consider a bipartite matching model [(C, 5, E), fi, Pol]. Assume that the structure 
(C, S, E, F) is stable, equivalently that {C U E U F) is strongly connected. Then the transition 
graph of the Markov chain of the bipartite matching model satisfies the property UTC. 

Proof of Theorem \4.3\ We are going to prove that the empty state can be reached starting from 
any state. This is a sufficient condition for property UTC to hold. The unique terminal strongly 
connected component is the set of states that can be reached from the empty state. 

We carry out the proof in the commutative case, but it works unchanged in the non-commutative 
case (the only information needed is the number of customer/server of each class). Consider a 
non-empty state (X, y), with X = (xc)cec and Y = {ys)s^s- It is sufhcient to prove that we can 
always reach a state {X' ,Y') such that \X'\ < \X\. 

If there exists (c, s) S C@ x S@ such that fi{c, s) > 0, then the proof is completed. Assume 
now that /i(C@ x S"®) — 0. Choose (c, s) S C@ x S*®. By assumption, there exists a path from 
s to c in {C U S, E U F). Let us denote it by (s = si, ci, S2, C2, . . . , S/t, = c). Assume that 
ci, . . . , Cfc_i ^ C@ and S2, . . . , Sk ^ S@. (If not, consider a subpath with this property.) Assume 
further that ci G Co- (If ci G C,, then (ci,S2) G E implies S2 G C@, so we can consider the 
subpath (s = S2, C2, . . . , Sfe, Cfc = c).) Since {si,Ci) G F and (c,;,Si+i) G by construction and 
since ci G Co, we get that ci, . . . , c^-i G Co and S2, ■ • ■ , Sfc & So- 




Figure 10: The path (s = si, ci, S2, C2, . . . , Sfe, Cfc = c) in {C U S,EU F)- 

By definition of the graph (C U S, E Li F), we have n{ci, Si) > for all i. Choose the sequence 
of arrivals (ci, si), . . . , (cfc, Sfe). Consider the effect of the arrival of (ci,si). Since si G S'@, it 
will be matched with a customer of C, (and not with ci, even if (ci, si) G E, since an admissible 
matching policy is always buffer-first, see j |2.2[ ). Since ci ^ C@, it will remain unmatched. Let 
(X(i),y(i)) be the new state. We have \X^^^\ = \X\- Also, in the new state, we have ci G c!,^\ 
which implies that S2 G S'@''. So we can repeat the argument inductively. After the arrivals of 
(ci, Si), . . . , (cfc_i, Sfe_i), we are in a state (^x'-''^^\Y^''^^'>) satisfying: 

Therefore, after the arrival of (cfc, Sfc), we end up in a state {X'-'^^Y'-''^) such that = |A:|-1. 

This completes the proof. □ 



15 



Example 4. Consider a bipartite matching model associated with the structure (C, S, E, F) where 
(C, S, E) is the NN graph of Figure [2| and where 



F= 1(1,1'), (2, 2'), (3, 3')} 



4.2 



the graph satisfies prop- 



The graph (CDS', EUF) is strongly connected. According to Theorem ' 
erty UTC. But it is not irreducible. Indeed, it is impossible to reach the state ((0, 1, 0); (0, 0, 1)) 
starting from the empty state. More generally, none of the states of the facet ({2}, {3'}) belong 
to the terminal strongly connected component. 

Below, we study the stability of bipartite matching models. Therefore, we always assume that 
the necessary conditions NCond are satisfied. So we get the property UTC for the Markov chain 



as a consequence of Theorem 4.3 



5 Models that are stable for all admissible policies 

Definition 5.1. Consider a bipartite graph (C, S, E) and an admissible matching policy Pol. The 
stability region is the set of values of fi for which the bipartite matching model [(C, S, E), /i, Pol] 
is stable. 

The stability region is included in the polyhedron defined by NCond. The stability region is 
maximal if it is equal to this polyhedron. 

Denote by ^ the set of facets. Define the following conditions on /i: 



scond: fic{c@{J')) + ^ls{s@{J')) > i-t^{Enc,{j^)xS,{j^)), yj^es-m (12) 



Let be a saturated facet, see Definition |2.4| Assume for instance that Co{J-) = 0. Then 
E n Co (J") X 5'o(J") = and C®(J") ^ C - C,{T). So pi implies: 



Since <S'@(J^) — S{C,{J^)), we recognize exactly ([5| for U = C,{F). Conversely, consider U C. C 
and the associated condition in NCOND: fic{U) < fis{S{U)). Choose a state with a strictly 
positive number of customers/servers for the classes U and S — S{U). Let be the corresponding 
facet. The facet J" is saturated: S@{T) = S{U), S,{T) = S - S{U), Soi^) = 0. Let us apply 
( 12 1 to the facet J", we get: 

Mc(C@(J-)) + M5(5(C/)) > 1 



To summarize, the subset of the inequalities (12) obtained by considering only the saturated 
facets gives precisely the inequalities NCond. 

We now show that the conditions SCond are sufficient stability conditions. 

Proposition 5.2. A bipartite model with probability fi satisfying SCOND is stable under any 
admissible matching policy. 

Proof. Consider the linear Lyapunov function: 

L{u, v) ~ \u\, {u, v) £ £, 

the number of unmatched customers (servers). Let (C/„, Vn)n be the Markov chain of the buffer- 
content. Let ^ be an arbitrary and fixed facet. Then for any {u,v) G we have (see 
Table [l]): 

E[L(i7„+i,14+i) I {Un,Va) = {u,v)]~Liu,v) = -Ai(C®(J-),5®(J-)) + Ai(Co(-F),5.(J-)) 
+ fi{C,{T),S,iJ'))+fi{C,{T),So{J'))+KCo{J') X So{J')nE-) 

= 1 - /ic(c@(J-)) - ns{s@{r)) - ti{Co{r) X So{:f) n e) . 
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The inequaUty (12) imphes directly that: 

E[L{Un+l,Vn+l) I {Un,Vn) = {u, v)] ~ L{u, v) < £ < 0. 



(13) 



By application of the Lyapunov-Foster Theorem, see for instance [J §5.1], we conclude that 
the model is stable. □ 







Co 


c. 




-1 








So 





or 1 


1 


s. 





1 


1 



Table 1: Variation of the linear Lyapunov function. 



Corollary 5.3. Consider a bipartite graph in which any non-zero facet is saturated. For any 
admissible matching policy, the stability region is maximal. 

The bipartite graph (C = {1, 2}, S = {1', 2'}, CxS- {(2, 2')}) is such that any non-zero facet 
is saturated. Therefore, its stability region is maximal for any admissible policy. The same is true 
for the "almost complete graphs" {C — {1, . . . ,k}, S ^ {1', . . . , fc'}, CxS— {{i, i'), Vi}). 

Example 5. Consider the NN graph from Figure [2] The graph has only one non-zero facet that 
is non-saturated, facet ({3}, {3'}). For any admissible policy, the stability region is at least the 
polyhedron SCond, Proposition |5.2[ which is defined by: 



NCOND, Aic(l)+/is(l') > 1-^^(2,2'). (14) 
Assume now /i = x Ms and fic — Ms- Set x — /ic(l) — Ms(l') a-nd y = l^c{'^) = Ms(2')- Then: 

NCOND : I ^ , SCond : 

2x + y > 1 



NCOND 

2a: + > 1 



In Figure 11 the light (yellow) region corresponds to SCOND, and the union of the light and 




Figure 11: NCOND and SCOND for the NN-graph with ^ = fic x and fic = Ms- 
dark (red) regions corresponds to NCond. 

Unfortunately, for some bipartite graphs, the polyhedron SCond is empty. This is illustrated 
by the following example. 



Example 6. Consider the NNN graph of Figure 12 The condition SCond for facet ({1}, {4'}) 
gives: 

Mc({3, 4}) + fisiil', 2'}) > 1 - ^(2, 3') (15) 

and for facet ({4},{1'}): 



Mc(l) + Ms(4') > 1 - m(2, 2') - m(2, 3') - ^(3, 3')- 



(16) 
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The inequality (15) is equivalent to: Aic(l) + Ms(4') < 1 — /^c(2) — /i({l, 3, 4}, 3'). Together with 
( 16 1 this gives: 



Mc(l) + M5(4') < 1 - Aic(2) - 3, 4}, 3') < 1 - ^(2, 2') - fi{2, 3') - ^(3, 3') < fic{l) + Ms(4'), 
which is impossible. 





1' 2' 3' 4' 1' 2' 3' 

Figure 12: NNN graph: facets ({1},{4'}) and ({4},{1'}). 




{11,(4'1 (1),(3'1 {21,(4'1 (21, ID 131.12'1 13),{r| {41,{3'( {41,(2'| W,[V] 

(ii}.i3\4M )( uii.m) ( (i,4uTD ( j2Kn^^ 



C(2.3.41,ll'l) C[3,4),(r,2'0 041,|r,2',3'} ) 

Figure 13: Facets for the NNN graph. Saturated facets are encircled (13 among 25 facets). 



6 Priorities and MS are not always stable 

Consider the NN bipartite graph of Figure [2] and Example [5] For this model, Proposition |5.2 
does not allow to decide if the stability region is maximal (see Figure 11). In Figure 14 we give 



simulation results for the average buffer size up to time n = 1000000 for the NN-graph with 
fi = X Ms, ^J'C = Ms, and MS policy. We can see that the average buffer size is growing rapidly 




Figure 14: Average buffer size for the NN-graph with ^ = /^c x A*s, Mc = Ms, and MS policy. 

near the 2x + y = 1 line. This does not necessarily imply unstability, as even for stable models 
we could have the mean stationary buffer size that is growing unboundedly as we approach the 
boundary of the stability region. In fact, we show below that for the PR and MS matching policies, 
the stability region is not maximal. 
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Proposition 6.1. Consider the NN model with either the MS policy or the PR policy given by: 



2 1 
2 1 

1 



id B = 



2 1 
2 1 

1 



For both policies, the stability region is not maximal. 



Proof. We carry out the proof for the PR poHcy. The idea of the proof is to play with two different 
Markov chains: the one describing the evolution of the buffer content, and an auxiliary one which 
mimicks the evolution of the customers/servers of type 2/2' on some of the facets. 

Consider an auxiliary Markov chain on Z with transition probabilities: 

X ^ X — 1 X ^ X X ^ X + 1 



X <0 
x^O 
x>0 



b-i 

C-l 



flo 
Co 



bi 

Cl 



Assume that a_i, ai, &_i, 6i, c_i, ci, are all different from 0. The chain is positive recurrent iff: 
[a_i < fli, Cl < c_i]. The stationary distribution is then equal to: 



b-i 



a-i 
bi 



7r(0) 
Tr{x) 

nix) = .(0)'^("^) 



bi 



Cl Vc_i/ 
b-i /a-i 



-1 



\x\ 



C-l - Cl- 

if X > 0, 
if a; < , 



Consider a NN-model with a probability ^ such that supp(/i) ~ C x S. Let {X, Y) = 
{X{n),Y{n))n be the Markov chain of the buffer content, where X{n) = {Xi{n), X2{n), X^in)) 
and F(n) = {Yi{n),Y2{n),Y3{n)). Assume wlog that {X,Y) is given under the form of a Stochastic 
Recursive Sequence, that is: 

(Xin + 1), r(n + 1)) = y(n), 0„) , 

where (0„)„ is an i.i.d. sequence of r.v.'s distributed according to fi, and $ is a deterministic 
function. 

Consider now the process {X2{n) — I2 ("•))?!■ This is not a Markov chain. However, if X2{n) + 
X^in) > 0, then "it becomes one". More precisely, if X2{n) + Xr^{n) > 0, then 

X2{n + 1) - Y2{n + 1) = $2(^2H - Y2{n), e„) , (17) 

where $2 is a deterministic function. This can be checked by direct inspection. Moreover, the 
transition kernel on Z defined by the recursion (17 1 is of the type of the above auxiliary chain 
with parameters: 



ai 
ao 

a-i 
bi 
bo 

b-i 

Cl 
Co 
C-l 



/^(l,l')+A*(l,3') + /i(2,l')+M(2,3'), 
/i(l,2')+A^(2,2') + /i(3,l')+M(3,3'), 

/i(2,l')+A*(2,3'), 



m(i,i') 

m(1,2') 
m(2,3'), 
m(1,3') 



Mi,3')- 

Ai(3,2'), 

M2,i')- 



m(2,2')+m(3,1') + A^(3,3'), 



M2,2')+/i(3,3'), 
= m(1,1')+A'(1,2')+M3,1')+m(3,2'). 
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Let us justify for instance the values of c_i,co,ci. We are in the case X2{n) + X3{n) > and 
X2{'n) — Y2{n) > which implies: 

Xi{n) = 0,^2 (n) > 0, = ¥2(71) - 0,1^3 H > , 

In Table [2] below, we show the effect of the different possible types of arrivals, restricting to the 
ones which may affect X2, i.e. when the customer class is 2 or the server class is 1' or 2'. To 
simplify, we have assumed in Table [2] that X^ > 0. For X3 = 0, the "Possible matchings" column 
would be affected, but not the "Selected matchings" and "AX2" columns. 



Arrival 


Possible matchings 


Selected matchings 


AX2 


(1,1') 


(1,3'), (2,1'), (3,1') 


(1,3'), (2,1') 


-1 


(1,2') 


(1,3'), (1,2'), (2,2') 


(1,3'), (2,2') 


-1 


(2,1') 


(2,1'), (3,1') 


(2,1') 





(2,2') 


(2,2') 


(2,2') 





(2,3') 








+1 


(3,2') 


(2,2') 


(2,2') 


-1 


(3,1') 


(2,1'), (3,1') 


(2,1') 


-1 



Table 2: Effect of arrivals. 



Let us comment on a couple of cases. If the arrival is of type (1,1'), then the selected matching 
is (2, 1') rather than (3, 1') due to the PR policy (i?2,i' > ^3,i')- H 1^^^ arrival is of type (1, 2'), the 
selected matching is (2,2') rather than (1,2') according to the buffer-first property of admissible 
policies, see §2.2[ The other cases are argued similarly. 



Let us introduce a new Markov chain (Wn)n on Z defined by: 

W^„+l=$2(l^n,^?„). 

(The process (Wn)n is different from the process {X2{n) — l2(n))„. The former is always defined 
according to the recursion pTl while the latter is defined according to ([T7| only for the n's such 



that X2{n) + X:}{n) > 0. The former is Markovian while the latter is not.) 

Condition ci < c_i becomes /ic(2) < /is(l') + ^s(2') and a_i < ai becomes ^s(2') < fj.ci^) + 
/ic(2). Both conditions follow from NCond. So the auxiliary chain {Wn)n is ergodic and its 
stationary distribution tt satisfies: 

^(0)=(l+ + Y\ 7r(Zl) -7r(0)— , n{Z*_)^TT{0) , 

V ai — a_i c_i — Cl/ c_i — Cl ai — a_i 

where Z^ = {1, 2, . . .} is the set of strictly positive integers and Zl — {—1, —2, . . .} is the set of 
strictly negative integers. From now on, we fix an initial condition Wq satisfying 

Wo r^TT, WoAL {en)n ■ 

Let us switch back to the Markov chain {X,Y). Set 

L{n) = X2in) + X^in), AL{n) = L{n + 1) - L{n) . 
If L{n) > 0, then we check by direct inspection that: 

AL{n) ^ ^{X2{n) -Y2{n),e„) , 
where is a deterministic function. We have in particular 

dpf 

a = E[AL(n) |L(n)>0,r2(n)>0] = Ai(3,2')+/i(3,3')-M(l,l')-M(2,l'), 

/3 =f E[AL{n)\L{n)>0,X2in)^Y2{n) = 0] = /i(2, 3') + ;u(3, 2') + ^^(3, 3') - /i(l, l'), 

dpf 

7 °" E[AL(n) |L(n)>0,X2(n)>0] = /i(2,3') + ^(3,3')-Ml,l')-/^(l,2'). 
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Let us turn again to the auxiliary chain {Wn)n- By performing the computation, we get 
E[^{Wn,9n)] - n{Z*_)a + 7r(0)/3 + ^(Z;)7 . 
The Ergodic Theorem for Markov Chains, see for instance il, §3.4], gives: 

^ n — 1 

hm - V -^{Wn, On) = TT{Z*_)a + n{0)l3 + 7r(Z* )7 a.s. 
n n ^ — ' 

1=0 

Assume that n{Z*_)a + 7r(0)/3 + tt{Z*^)j > 0. Then we have: 

n-l 

hm V^'(Ty„,6'„) = +00 a.s. (18) 

i=Q 

Therefore, for each £ > 0, there exists i^r^ > such that 

ra-l 

p(minV^'(W„,0„) > -if,) > 1-e. (19) 

\n>l ^ — ^ / 

~ 4=0 

Let us switch back to the Markov chain {X{n), Y{n))n- Choose the initial condition {X{Q), 1^(0)) 
such that 

X2{0)~Y2iO)^Wo min(X3(0),r3(0)) =ife, (20) 

where ifg is defined in ( 19 1. By construction, on the event A = {min„>i J27=o ^i^n, On) > — i^e}, 
we have 

Vn, L{n)>0, X2{n) - ¥2(11) = Wn . 
So, on the event A, we have 

n — 1 n— 1 

L{n) ^K, + J2 ^L{i) ^K, + Y, d,) +00 . (21) 

i=0 i=0 

We conclude that the Markov chain {X, Y) of the NN-model is transient. 

We now show that the stability region is not maximal, by giving an example such that 7r(Z*L)a+ 
7r(0)/3 + 7r(Z^)7 > 0. Consider fic = (1/3, 2/5, 4/15), fis = /^c, and fi = fi^x Thus conditions 
NCOND are satisfied. However, we have: 





11 




34 




8 


ai 


" 25' 




75' 


a-i = 


75' 




6 




13 




6 


hi 


" 25' 


60 = 


25' 


6-1 = 


25' 




8 




34 




11 


Cl 


~ 75' 


Co = 


75' 


C-l = 


25' 



and 

This gives a = -1/15, P = 13/75, 7 = -1/15, and, 

29 
915 

This completes the proof. 

Consider now the MS policy. Set L{n) = min(X3(ri) — X2(n), Y-i{n) — Y2{n)). The initial 
distribution can be taken such that 



TT{Z*_)a + 7r(0)/3 + 7r(Z;)7 = — > . 



^2(0)- 1^2(0) -TT, 



|X3(0)-X2(0) =Xe ifX2(0)>0 

[r3(0) - ^2(0) = Ke otherwise 
Modulo these modifications, the proof carries over unchanged. □ 
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7 ML is always stable 

In this section, we show that the ML pohcy has a maximal stabihty region. 

The idea of the proof is as follows. Consider the quadratic Lyapunov function: 

L{x,v)^Y.^^c + Y.y^s, (a^,y)e£. (22) 

Observe that the ML policy minimizes the value of this Lyapunov function at each step. We 
introduce an alternate policy that depends on the arrival distribution /z. For this policy, we 
manage to prove that the quadratic Lyapunov function has a negative drift outside a finite region. 

Theorem 7.1. For any bipartite graph, the ML policy has a maximal stability region. 

Proof. We introduce an alternate matching policy. This policy is admissible, corresponds to a 
commutative state space, but does not belong to the policies listed in §2.2| It is a random policy 
and its specificity is to be facet dependent. 

Let us describe the alternate policy on a non-empty facet J^. Set C, — C, ( J^) , S, = S, ( J^) , C@ = 
C(S){T), etc. To describe the matching policy, the only thing we have to describe is the way to 
match an arriving customer of class c € C@, resp. server of class s £ S@. Let us concentrate first 
on a server of class s £ S*® . 

From NCOND: 

We build a directed graph as in ^ but restricted to the nodes in C, and S@. Formally, 

A(f = (C. U S*® U {^, /}, {EnCx 5®} U {{i, c), c G C,} U {(s, /), s e 5®}) . (23) 

Endow the arcs of i? H C, x 5® with infinite capacity, an arc of type (i,c) with capacity ^,c{c), 
and an arc of type (s, /) with capacity usis). 



As in Lemma [3.3[ NCond implies that the minimal cut of A/jp has capacity /ic(C,). Any 
maximal flow T is such that: Vc e C,, T{i,c) = /ic(c), Vs £ S"®, T{s,f) < /is(s). Let us prove 
that there exists a maximal flow T such that: 

Vc e C,, T{i,c) = ncic) 
Vse5®, T{sJ)<ns{s)- 

Define Jls on S@ by Jls{s) — l^s{s) ~ V- Here 77 > is chosen to be small enough so that: 
V[/ C C,,nc{U) < Jls{S{U)), yV C S@,JlsiV) < nc{C{V)). This is possible since NCond 
are open conditions. Consider the same network as above but with the capacities Jls{s) on the 
arcs (s,/). The minimal cut still has capacity ^i.c{C»). A maximal flow T is such that: Vc G 
C,, T{i,c) = lJ'c{c)^ Vs € S(Sn T{s,f) < Jls{s) < fJ-s{s)- Clearly, T is also a flow for the original 
network. 

The server s £ 5'® is matched to c € C, n C(s) randomly, independently of the past, with 
probability: 

1 



Ais(s) 



\c,nc{s) 



Let us check that this deflnes indeed a probability: 



cec.nc(s) ' cec,nc(s) 



1 



fisis)-TisJ)+T{sJ) ] 



22 



For ceC, s £ S{c), set e^c = {f^s{s) ~ T{s,f))l\C, n C{s)\. For c G C., set = Eses(c) ^^c- 
We have £c > 0. Observe that: 

VceC, 5] A'5(5)P,-^ = Mc(c)+ec. (24) 
ses(c) 



SymmetricaUy, we define the directed graph of type ( p3| but on the nodes C@ and S, . We build 
a maximal flow on this new graph as above, and based on this flow, we define the probability 
that a customer c € C@ is matched to a server s e S*. n S'(c). For s € S',, we define Scs, c e C(s), 
and Es accordingly. We have Eg > 0. 

Let {X{n),Y{n))n be the Markov chain of the buffer-content of the model. Assume that 
{X{n), Y{n)) ^ {x, y) e F and let c e C,. We have: 

(i) X(n+1), = X(n),-liff: 

- the arriving customer is not of class c; 

- the arriving server is of class s € <S'(c); 

- the arriving server is matched with c (probability -P^). 
This case happens with probability ac = X]ses(c) /^(C* — c, s)P^. 

(ii) X(n+ l)c = X(n)c + 1 iff: 

- the arriving customer is of class c; 

- the arriving server is not matched with c. This may occur in two possible ways: either 
the arriving server is of class s ^ S{c), or the arriving server is of class s e S{c) but is 
not matched with c (probabihty 1 — P^)- 

This case happens with probability /3c = /i(c, 5* — 5(c)) + J2ses{c) /^("^' ^)(1 ^ Psc)- 

(iii) Otherwise, X{n + l)c — X{n)c- 
Using (p4|), we get: 



= /xc(c)+ec- E A'(c,s)-Ps^ 

= /i(c,5-5(c)) + /i(c,5(c))+ec- E Ai(c,s)P,-^ 

ses(c) 

= ^(c,5-5(c))+ E /i(c,s)(l-P,-^)+£,. 

s<£S{c) 

Thus, Qfc = /3c + £c- Observe that (3c < Mc(c)- We get, for {x, y) £ T and c G C,: 

E[X(n + l)^-X(n)2 I (X(n),r(n)) - {x,y)\ = /3c(2.Tc + 1) - ac(2a:c - 1) 

= 2/?c - £c(2a;c - 1) 



Let L be the quadratic Lyapunov function (22 1. Define AL{n) = L{Xn+i,Yn+i) — L(X„,y„). 
Set e = nrin^gc^us. > 0. Then (the first term in the sum takes care of the vertices in C — C, 
and S ~ S,): 

E[AL(n) I (X„, y„) = {x, y)]<2+Y, C^f^cic) ~ e^Xc) + ^ (2Ms(s) - e^y.) 

cgc. ses. 



< 2 + 2^c(a) + 2Ais(5.) - £ ( E + E ^ 



< 6 - 2 £ E 
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Fix S > 0. If J2cec '^c > (6 + S)/2e, then E[AL(n)] < —5. There are finitely many facets, so there 
is a finite set A C £ such that 



y{x,y)^A, E[AL{n)]<-5. 



(25) 



By the Lyapunov-Foster's Theorem, see for instance \Tl §5.1], the alternate matching policy is 
stable. 

Since the ML matching policy minimizes the value of the quadratic Lyapunov function, we 
have a fortiori that (25) holds for it. Therefore, the ML policy is also stable. □ 



Conclusion. Many open questions remain. First, we do not know if the stability region is always 
maximal for the FIFO and Random policies. Numerical experiments seem to indicate that it is 
indeed the case. Second, for the MS and priority policies, we know that the stability region is not 
always maximal, but we do not know how to compute it. Last, we would like to obtain sufficient 
conditions for stability, valid for all admissible policies, and which are better than the ones of ^ 
The program used to carry out the numerical experiments is available on request from Ana 
Busic. 
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